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Abstract 

The interferometric analysis of meson correlations provides a measure of the average phase 
space density of the mesons in the final state. This quantity is a useful indicator of the 
statistical properties of the system, and it can be extracted with a minimum of model 
assumptions. Values obtained from recent measurements are consistent with the thermal 
value, but do not rule out superradiance effects. 



It would be interesting to know the average phase space density of the pions produced 
in ultrarelativistic heavy ion collisions. In the final state, the local phase space density is 
frozen (Liouville's theorem) and it gives a measure of the dynamics in the prior interacting 
region. If the system could be described by a local statistical equilibrium, the distribution 
function / would have the Bose-Einstein form, frip^f) = l/(exp(u(r) ■ p/T) — 1). For 
massless particles the average of this quantity is a pure number, 

r d^pd^rf^ C(2) 
J d-^pd-^rfr C(3) 

For massive bosons, the number is lower; for example, for pions at a chemical freezeout 
temperature of 200 MeV the average is 

(/)r.,T~0.15. (1) 

If the experimental (/) were close to this, it would be welcome evidence for the existence 
of a local equilibrium. If (/) came out much larger, it would lend considerable support to 
the idea of superradiant pion states[l], which have been an object of renewed interest [2, 3]. 
On the other hand, if (/) came out much smaller than thermal, it would point to entropy- 
generating processes such as the slow decay of heavy resonances or quark-gluon droplets. 
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In this letter I want to point out that the measured two-particle correlations yield 
direct information about the phase space density of the mesons, when interpreted according 
to Pratt's interferometric formula[4]. To see this, we begin from Pratt's formula written 
as 
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where g is the source function for the mesons. We next convert the source function g to 
an equivalent source at a common time Iq by the replacement 

g{r, t,p)^ 6{t -to) f° dt'g{f- v{t' - to), t', p) = S{t - to){27r)^f{r,p) 



where v is the velocity associated with the momentum vector p. This replacement does 
not affect the correlation function if w • (pi — p^) = {Ei — E2). The condition is satisfied for 
small momentum difi'erences and for longitudinal motion of extreme relativistic particles, 
and seems rather safe for the present application. The coefficient of the 5-function is then 
the phase space density at to, extrapolating the positions of the final state mesons to 
that time. We next integrate over the momentum difference d^q, which produces a delta 
function 5^{ri — r2) to eliminate one of the spatial integrals. The result is 

J J d'^xid'^X2g{xi,p)g{x2,p) cosq- {xi- X2) ^ {27t)^ J d^rf{r,p). 

Finally we integrate over d^p/ {27r)^ and normalize to the number of particles to obtain the 
phase space average. 
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For heavy ion collisions it is more useful to make the average over a fixed rapidity interval, 
because the system evolves to produce a spatial separation between particles of different 
rapidities. The different rapidity groups equilibrate independently. The formula for a small 
rapidity interval reads 



d'^Ptidyid'^Pt2dy2 p+q,p-q d'^ptdy p+qd'^ptdy 
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Eq. (3) and (4) just require integral properties of the correlation function, so they should 
be less dependent on the accuracy of the momentum measurements than other observables. 

The integral should undoubtedly be evaluated directly from the experimental data, 
but for an orientation I shall try to evaluate it from published Na44 parameterized distri- 
butions [5]. The Na35 experiment [6] obtained similar information, but did not quote the 
entire parameterization. One of the common parameterizations for these correlations is 
Gaussian with parameters A and source sizes Rl,Rs and Ro, 

.2 , ,2 , = 1 + Aexp + q^Rl + qlRl) —— . 

d'^Ptidyid^Pt2dy2 \ 2 J d-^ptdy pi d^ptdy P2 

The single-particle transverse momentum spectra can be parameterized by an exponential 
function, 

d^n^^^ dn exp{—pf /Tt) 
d^ptdy " Ity 

With this parameterization, the average phase space density is given by 
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RhRaRoTf 

Ref. [5] quotes the following numbers for S + Pb— > 7r"''+X at midrapidity: A ^ 0.4, 
Rt ~ 6.0 fm, Rl ~ 6.0 fm and dn/dy fa 40. From their Fig. 7 can be deduced Tt ^ 187 
MeV/c. Taking Rg = R^ — R^ and inserting these numbers in eq. (3), I obtain 

(/)d2/~ 0.07 -0.16. 

The range is obtained from the quoted experimental errors combined quadratically. I note 
that the NA35 experiment found a somewhat smaller source size; the lack of agreement 
between experiments is a caution not to draw strong conclusions from the present data. 

One should also be reminded that the assumptions going into the interferometric 
formula, eq. (2), may not be well satisfied. The most critical assumption is that the one- 
particle distribution function is not affected by the Bose symmetrization, which is only 
satisfied for low phase space densities. 

Given these caveats, what does one conclude? The extracted phase space density is 
consistent with local statistical equilibrium for a chemical freezeout temperature in the 
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100-200 MeV range. The analysis includes pions from long-lived resonance decays such as 
a; ^ Stt and 77 Sir. If these could be subtracted out, the phase space density would be 
somewhat higher. Thus there may be an excess of pions produced in the source, and the 
possibility of coherent pion effects should not be ruled out. 

The extracted phase space density appears high enough to make unlikely that long- 
lived intermediates such as droplets of quark-gluon plasma are produced abundantly. This 
is consistent with the present theoretical expectation of no strong first-order phase transi- 
tion in the quark-gluon plasma [7]. 
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